Gauge fluctuations and transition temperature for superconducting wires 
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We consider the Ginzburg-Landau model, confined in an infinitely long rectangular wire of cross- 
section Li x Z/2. Our approach is based on the Gaussian effective potential in the transverse unitarity 
gauge, which allows to treat gauge contributions in a compact form. The contributions from the 
scalar self-interaction and from the gauge fluctuations are clearly identified. Using techniques from 
dimensional and zefa-function regularizations, modified by the confinement conditions, we investi- 
gate the critical temperature for a wire of transverse dimensions Li, L^. Taking the mass term in 
O ' the form m\ = a(T /Tq — 1), where To is the bulk transition temperature, we obtain equations for the 

»' | , critical temperature as a function of the L[s and of To, and determine the limiting sizes sustaining 

f> i 1 the transition. A qualitative comparison with some experimental observations is done. 
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£H , I. INTRODUCTION 

o 

In the present work we discuss the critical behavior of the Ginzburg-Landau model compactified in spatial dimensions 
and we implement the gauge contributions using the Gaussian effective potential formalism. Spontaneous symmetry 
breaking is obtained by taking the bare mass coefficient in the Hamiltonian parametrized as uiq = a(T/T Q — 1), with 
a > and the parameter T varying in an interval containing the bulk transition temperature T . With this choice, 

■ considering the system confined in an infinitely long rectangular cylinder with cross-section area A = L\ x i 2 , we 
obtain the Ginzburg-Landau model describing phase transitions in samples of a material in the form a wire. This 

OO 1 generalizes to wires recent results obtained for materials in the form of films [1]. We investigate the behavior of 
the system taking into account gauge fluctuations, which means that charged transitions are included in our work. 
We are particularly interested in the problem of how the critical behavior depends on the dimensions L\ 1 L2- This 
study is done by means of the Gaussian Effective Potential (GEP) as developed in Refs. [2-5], together with a spatial 

■ compactification mechanism introduced in recent publications [8,9]. 

A central ingredient in our approach relies on the topological nature of the Matsubara imaginary-time formalism. 
I ■ In order to perform the calculation of the partition function in a quantum field theory, the Matsubara prescription 
results to be equivalent to a path-integral calculated on R ^ 1 x S 1 , where S 1 is a circle of circumference (5 = 1/T. 
This result was demonstrated at the one-loop level in Ref. [10] and has been assumed to be valid for higher orders [11]. 
Relying on the fact that for Euclidian field theories inverse imaginary time and spatial coordinates are on the same 
footing, such a topological result has been generalized to treat different physical situations in which fields are confined 
in higher purely spatial dimensions, considering the Matsubara mechanism on a R D ~ d x 5 11 x 5 l2 ... x S ld topology, 
describing space confinement in a <i-dimensional subspace [8,9,12]. As a consequence the Matsubara formalism can 



in 
o 



o 
o 



be thought, in a generalized way, as a mechanism to deal with spatial constraints in a field theory model. In this 
situation, for consistency, the fields fulfill periodic (antiperiodic) boundary conditions for bosons (fermions). In any 
case we infer from the above discussion that we are justified to consider in this paper the Matsubara mechanism as a 
path-integral formalism on R D ~ 2 x S 1 x S 1 to deal simultaneously with spatial constraints in a subspace of dimension 
d = 2 (the transverse dimensions of a wire having a rectangular cross-section). These ideas have been applied in 
different physical situations [8,12,9,13]. 

We begin by briefly presenting the study of the U(l) Scalar Electrodynamics in the transverse unitarity gauge, 
along the lines developed in refs. [1,6], and afterwards presenting a general procedure to treat a massive field theory in 
a D-dimcnsional Euclidean space, compactified in a c£-dimensional subspace, with d < D. This permits to extend to 
an arbitrary subspace some results in the literature for the behavior of field theories in presence of spatial boundaries 
[9,8,7]. After describing the general formalism, we focus on the particularly interesting case of d = 2. We investigate 
the behavior of the system taking into account gauge fluctuations, which means that charged transitions are included 
in our work. We are particularly interested in the problem of how the critical behavior depends on the dimensions 



Lj. This study is done by means of the Gaussian Effective Potential (GEP) as developed in Refs. [2-5], together with 
the spatial compactification mechanism mentioned above, introduced in recent publications [8,9]. 



II. THE GAUSSIAN EFFECTIVE POTENTIAL FOR THE GINZBURG-LANDAU MODEL 



We start from the Hamiltonian density of the GL model in Euclidean (i-dimensional space (unless explicitly stated 
we work in Natural Units (NU), h — c = k = 1) written in the form [14], 



-F, V F^ + \\{d»- ^)vl/| 2 + im^| 2 + A(|*| 2 ) 2 , 



(2.1) 



where ^ is a complex field, and m is the bare mass. The components of the transverse magnetic field, = 
d^Av — d v A u {n, v = 1, d) are related to the d-dimensional potential vector by the well known equation, 



1 



FauF^ = IV x A| 2 . 



(2.2) 



In order to obtain only physical degrees of freedom, we can introduce two real fields instead of the complex field "f, 
assuming a transverse unitarity gauge. We can define the field in terms of two real fields, as W = 0e* 7 , together with the 
gauge transformation A — > A— l/e'V'y. The unitarity gauge makes the original transverse field to acquire a longitudinal 
component proportional to V7. Then the original functional integration over ^ and ty* in the generating functional 
of correlation functions, becomes an integration over <f>, At and A^. The longitudinal component of the vector 
potential can be integrated out, leading to the generating functional (up to constant terms), 



Z[j] = J D(t>DA T ]exp - J d d xH + j d d x j<j> 



where the Hamiltonian is 



H 



2 (W 



-e 2 4> 2 A 2 ' ^ 2 ' 1 ^ » : 

2 



-(vxAr+-(v-Ar. 



(2.3) 



(2.4) 



We have introduced above a gauge fixing term, the limit e — > being taken later on after the calculations have been 
done. In Eq.(2.4) and in what follows, unless explicitly stated, A stands for the transverse gauge field. 

The Gaussian effective potential can be obtained from Eq.(2.4), performing a shift in the scalar field in the form 
4> = <j> + <p, which allows to write the Hamiltonian in the form 



H — Ho + Hint, 

with Ho being the free part and Hint the interaction part, given respectively by 



Ho — 



i(V0) 2 + i^ 2 



+ 



i(VxA) 2 + iA 2 V4^ + l(V-A) 2 



and 
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H in t = ]T vj n + - (eV - A 2 ) A U A» + -e 2 4>A u A»<p + -e 2 A^A^ 2 



n=0 



where 



(2.5) 



(2.6) 



(2.7) 



vi = rriQip + 4A<p 3 , 
V2 = \m 2 oV 2 + 6A<p z 



-n 2 



v 3 = 4X(p, 
v 4 = A. 



(2.8) 
(2.9) 
(2.10) 

(2.11) 
(2.12) 



It is clear from Eqs. (2.5, (2.6) and (2.7), that H describes two interacting fields, a real scalar field <f> of mass and 
a real vector gauge field A of mass A. 

The effective potential, which is defined by 



V, 



eff[ 



1 

V 



-InZ + J d d xjip + J d d xJ v A v 



(2.13) 



where V is the total volume, can be obtained at first order from standard methods from perturbation theory. One 
can find, from Eqs. (3.1), (2. 6) and (2.7), 



V e ff[<p] = J?(fi) + 2/f(A) + X -m 2 <p 2 + \cp 4 + 1 [m 2 - tl 2 + I2\p 2 + 6XI d (n)] I d (n) 
+ [e 2 (^ 2 + / d (f7))-A 2 ]/ d (A), 



where, 



and 



Io(M) = J 



d d k 



1 



(27r) d k 2 + M 2 ' 



I d (M) 



d d k 
{2ir) d 



ln(fc 2 + M 2 ), 



(2.14) 



(2.15) 



(2.16) 



with k = (fci, kd) being the d-dimensional momentum. 

The Gaussian effective potential is derived by the requirement that must be stationary under variations 

of the masses A and fl This means that values and A for the masses fi and A should be found such that the 
conditions, 



an 2 



= o, 



(2.17) 



dv e 



eff 



dA 2 

be simultaneously satisfied. These conditions generate the gap equations, 



= 0, (2.18) 

A 2 =A 2 



Q = ml + 12\tp 2 + 12A^(Q) + 2e 2 /^(A), (2.19) 
A = e 2 ^ 2 + e 2 I d (U). (2.20) 

Replacing £1 and A in Eq.(2.14) by the solutions f2 and A, of Eqs. (2. 19) and (2.20) we obtain for the GEP the formal 
expression, 

V ef f[<p] = lf(Q) + 2l((K) + \m 2 ^ 2 + \^ 3\[I d (n)} 2 e 2 I d (n)I d (K). (2.21) 

Notice that Eqs. (2. 19) and (2.20) are a pair of very involved coupled equations, and no analytical solution for them 
has been found, they can be solved only by numerical methods. Later on we will see that this difficulty, in the limit 
of criticality can be surmounted. 

Next we intend to write an expression for the Gaussian mass, m, obtained in our case from the standard prescription, 
as the second derivative of the Gaussian effective potential for ip = 0. To calculate the second derivative of V e ff with 

2 2 

respect to ip, we remark from Eqs. (2. 19) and (2.20) that f2 and A also depend on ip according to the relations 

dT i 2 24A^- e 2 ^ l( A)fl 
dip 1 + 6X1^(0) ' 

2 2 

n 9 1 9tW ,-A-sdQ, , n 

w =2eip ~2 eI - i{A W (2 - 23) 



where 

Replacing Eq.(2.23) in (2.22) we get, 



dfl _ [24A - 2e 4 /f 1 (A)] tp 
dtp ~ 1 + [6A - ±e 4 /f ^O)] /^(H) ' 

and the second derivative of the GEP with respect to tp is given by, 

' e// -ml + UXtp 2 + 12A/ d (n) + 2e 2 / d (A) + 2eV-i(A) 

[6A + ±e 4 /f ^A)] [24A - 2e 4 Jf 1 (A)] <^ 2 
1+ [6A- ±e 4 J_i(fi)] /^(fl) ' 



(2.25) 



# 2 



(2.26) 



Thus we have the formula for the Gaussian mass, 

= m 2 + 12A/ d (n o ) + 2e 2 / d (A ), (2.27) 



-2 d 2 V eff 

m = 



dtp 2 



ip=0 



where fl and A arc respectively solutions for and A at tp = 0, explicitly, 

sf = m 2 + 12A/ d (n ) + 2e 2 / d (A ), (2.28) 
A^=e 2 / d (f7 ). (2.29) 

Therefore, from Eqs.(2.27) and (2.28) we get simply, 

to 2 = Hq. (2.30) 



Hence, we see from the gap equation (2.27) that the Gaussian mass obeys a generalized "Gaussian" Dyson- Schwinger 
equation, 

to 2 = m 2 + 12A/ d (m) + 2e 2 / d ^Je 2 lg(m)j . (2.31) 
This expression will be used later to describe the system in the neighbourhood of criticality. 

III. CRITICAL BEHAVIOR OF THE GINZBURG-LANDAU MODEL COMPACTIFIED ON TWO 

SPATIAL DIMENSIONS 

We will work in the approximation of neglecting boundary corrections to the coupling constant. A precise definition 
of the boundary- modified mass parameter will be given later for the situation of D = 3 with d = 2, corresponding 
to a wire of rectangular section L\ x L 2 . Wc also consider the limiting situation of one of the transverse dimensions, 
L 2 — ► oo, which corresponds to a film of thickness L\. We use Cartesian coordinates r = [x\, ...,xa,z), where z is a 
(D — d)-dimensional vector, with corresponding momentum k = (fci, kd, q), q being a (D — c?)-dimcnsional vector 
in momentum space. Then the generating functional of correlation functions has the form, 

Z = J Vp^Vpexp I- J d d r j d D - d zH{p,\7p)\ , (3.1) 



where L = (Li, Ld), and we are allowed to introduce a generalized Matsubara prescription, performing the following 
multiple replacements (compactification of a d-dimensional subspace), 



+ oc 



Li ' Li 

rii — — oo 



, i = 1,2, ...,<i. 



(3.2) 



A simpler situation is the system confined simultaneously between two parallel planes a distance L\ apart from one 
another normal to the Xi-axis and two other parallel planes, normal to the a^-axis separated by a distance L 2 (a 
"wire" of rectangular section). We emphasize however, that here we are considering an Euclidean field theory in D 
purely spatial dimensions, we are not working in the framework of finite temperature field theory. Temperature is 
introduced in the mass term of the Hamiltonian by means of the usual Ginzburg-Landau prescription. 

For our proposes we only need the calculation of the integral given in equation (2.15) in the situation of confinement 
of the present section. With the prescription (3.2), the expression corresponding to equation (2.15) for the confined 
system can be written in the form, 



I D (M) 



1 



4ir 2 L 1 L 2 ..L d 



E 



ni,..,nd— — oo 



jD-d r 



(a 2 m\ + a\n\ + .. + a\n\ + c 2 + q 2 ) 



(3.3) 



In the following, to deal with dimensionless quantities in the regularization procedures, we introduce the parameters 



m t 

C= 2V ai = T i 



hi 
2^ 



(3.4) 



and using a well-known dimensional regularization formula [15] to perform the integration over the (D — d) non- 
compactified momentum variables, we have, 



iD-d 



Q 



n^T(s-^) 



1 



{a\n\ + a\n\ + .. + a 2 n 2 + c 2 + q 2 Y 
with which Eq.(3.3) is written, 



r( s ) 



(ami + a ^ n \ + ..+ a 2 d n 2 d + c 2 ) s 



(3.5) 



E 



4 7 r 2 L 1 i 2 ...L d (am 2 + a 2 n 2 + ■■ + a 2 d n 2 + c 2 )*-^ 



(3.6) 



The sum in Eq.(3.6) is one of the Epstein-Hurwitz zeta-functions, defined by, 

oo 

A c d {v;a\,...,a d )= H 1" a d n d + ( ?Y V 



(3.7) 



ni,...,n ti =-oo 



where v = s — . Next we can proceed generalizing to several dimensions the mode-sum regularization prescription 
described in Ref. [16]. This generalization has been done in [8] and we briefly describe here its main steps. ^From the 
identity, 



1 



1 



a- r(iz) J 

and using the representation for Bessel functions of the third kind, K v , 

2(a/b)* K v (2\fab) = / dx a^e"^" 6 *, 
Jo 

we obtain after some rather long but straightforward manipulations [8], 



(3.8) 



(3.9) 



A° d (v;ai, ...,a d ) = 



2»-i+ 1 ir 2 »-i 
ax---a d T(v) 



T-%- Y T\v 



(2ncf-^ + 2j2Y: 

i=\ nt = l 



2-Kca, 



K„ 



( 27rcn, 
V a i 



+ • • • + 2 d 



E 

ni,...,n d - 



=i V 



27rc V a 2 



K„ 



271-Ci / % 



n 2 d 



(3.10) 



In our case s = 1 and for wires d = 2, and Eq.(3.10) becomes, 



A 2 (2- T ,L 1> £ a )_ r(2 ^ #) 



2 -r ( :i - -)M 



D-2 



' t1 \ 00 / r2 x ^~ ~ 



ni — 1 

oo 



no — 1 



M 



K 1 _D(ML 2 n 2 )+ 



+2 2 E 

ni ,n 2 = l 



2^2 



1 — — 

L 1 



K x _d M\L{n{ + L{n{ 



(3.11) 



Thus, inserting the equation (3.11) into equation (3.6) and making some algebraic manipulations, we obtain 



ji\ = 1 

CO 



mi 



M 

where K v are the Bessel function of third kind. 



(3.12) 



IV. CRITICAL BEHAVIOR FOR WIRES 



In our case, we are analyzing a superconducting type-I wire, which obeys the following conditions [14], 

£(T) = (m)- 1 » A(T) » L, 



(4.1) 



where L is the smallest one of the linear dimensions in the rectangular cross-section area of the wire A — L\ x L 2l 
and £(T) and A(T) are respectively the critical correlation length and the critical penetration depth defined by, 



£C0 = 



t|V2 



; A(T) 



Ao 

t|V2 



t = 



T-T c 



(4.2) 



where T c is the transition temperature, £o and Ao the intrinsic coherence length and London penetration depth, 
respectively. In the case where one of the linear dimensions is much larger than the other, L\ > L 2 or L 2 > L\, we 
should retrieve a film-like behavior, and L in Eq.(4.1) would be the thickness of the film. After that, we can take 
M = to in Eq.(3.12) and restrict ourselves to the neighborhood of criticality, that is to the region defined by to m 0. 
Then we can use the asymptotic formula, 



#M(*) = k(M)(?) H ; 



which allows, near criticality, to write equation (3.12) in the form 



7 ( f(TO«0) = 2-* 7 r-*r(|-l) 



E 



2- 



or, 



/ D (m«0) 



_ ni =i (ii«i) 2 

r(#-i 



IPI+E 



2- 



i {L 2 n 2 ) 2 1 



2 E 



2T-1 



2^* 



if 



D-2 



-D-2 



C(£> - 2) + 2£ 2 ( 



L>-2 



; Li,L 2 ) 



where £7 2 (^pi Li,L 2 ) is the generalized 2-dimensional Epstein zeta-function defined by [27] 



E 2 (^-M,L 2 \= £ [Llnl + Ljni]-^ . 

^ ' ni , ri2 — 1 



(4.3) 



(4.4) 



(4.5) 



(4.6) 



In an analogous way as it is done for the Riemann zeia-function, one can also construct analytical continuations 
(and recurrence relations) for the multidimensional Epstcin-Hurwitz zeia-functions which permit to write them in 
terms of Kelvin and Riemann zeta functions. To start one considers the analytical continuation of the Epstein-Hurwitz 
zeia-function given by [16] 



1 r~ 

(« +P ) =-oP 



n=l 



n=l 



(4.7) 



Using this relation to perform one of the sums in (4.6) leads immediately to the question of which sum is firstly 
evaluated. Whatever the sum one chooses to perform firstly, the manifest L\ <-» L2 symmetry of Eq. (4.6) is lost; in 
order to preserve this symmetry, we adopt here a symmetrized summation. Generalizing the prescription introduced 
in [8], we consider the multidimensional Epstein function defined as the symmetrized summation 

00 00 

E d (v, L U L d ) = i ]T ]T • • • £ [alnl + ■■■ + a\n d \ ~" . (4-8) 

a ni = l 11,1 = 1 

where Oi — cr(Li), with a running in the set of all permutations of the parameters L\, ...,Ld, and the summations 
over m, ...,n d being taken in the given order. Applying (4.7) to perform the sum over n d , one gets 

1 d ~ 
Ed {v,L\, ...,L d ) = - — y^-Ed-i (v; .-.,Li, 

i—l 

+^ r H)g^- K-S-) + 5rf>»44") ■ (4 - 9) 

where the hat over the parameter Li in the functions Ed-\ means that it is excluded from the set {Li, ...,Ld} (the 
others being the d — 1 parameters of E d -i), and 

W d (rj; Ll ,...,L d ) = J2^ £ ( ™^ ) K J^y/{--. +£& + •••)) , (4.10) 

i=l z n 1 ,...,n d =i\L i ^/(--- + L 2 n 2 + ■■■)/ V 1 / 

with (• • • + L 2 n 2 + • • •) representing the sum Y^j=i E"]n 2 — E\ n l- ^ n particular, noticing that E\ (v; Lj) = LJ 2u ((2iy), 
one finds 

E, (^;Ll Lt)= -\ (' + ') ((D 2) 



4 Vif- 2 if" 2 



Inserting (4.11) into (4.5) we get, 

-»> ^ "i* 2 (ip + zH « D - 2 » + (i^ + sip) « D - 3 » 

r n D-l/2 ' 

where, 

W 2 (0;Li,L 2 ) = £ (j-KoUn^n^+^KoUnj-n^Y (4.13) 

rai, 712=1 1 2 2 

For D = 3, the first and second terms between the brackets in Eq. (4.12) are divergent due to the poles of the (^-function 
and T-function, respectively. However, the whole expression in (4.12) is not divergent, in fact these divergences cancel 
between themselves leaving a finite result for ^(m w 0). This can be easily seen if we consider the following properties 
of the zeta-function, 



C(*) = r^r(^V-i C (i-*) ; lim 



C(*)- 



z - 1 



= 7, 



(4.14) 



where 7 w 0.5772 is the Euler-Mascheroni constant. 

The quantity M^(0; L\,L 2 ), appearing in Eq.(4.12), involves complicated double sums; in particular it is very 
difficult to be analytically handled for L\ 7^ L 2 , what means that it is not possible to take analytically limits such as 
L 2 — > 00 for finite L\. Such a limit (which corresponds to a film of thickness L\) will be considered numerically. 
Using Eq. (4.14) in Eq. (4.12) we obtain for D = 3, 



7 



l§(m*0) « ^ I — + — ) + -W 2 (0;LiL 2 ) 



1 



(4.15) 



and Eq.(3.12) becomes, 



^o 3 (Ao) 



+ 



1 



£H)Ao 

23^3/2 1 2 1 /2 7r 3/2 
A 



E 



An 



1/2 



A"i /2 (AoLini) + 

«2 = 1 



+ E 



1 \yJL\n\ + L\n. 



Lint 

1/2 

K 1/2 (A ^J Ljnl + L 2 2 n 2 2 ) 



An 



1/2 



i^l/2(A 0j L2«2) + 



(4.16) 



Using the formula, 



Ky 2 {z) = \l^e- 



(4.17) 



and remembering the notation, Ao = \J e 2 / D (m) we obtain, after some rearrangements, 



I*(y/e*I$(m)) = - 



47T 



27Til 



2ttL 9 



+ 



1 ~ e - ev /^pK«0)Vi?n?+L|n| 

+^ E 



- n u n 2 = l V L i n i + L ln'l 

Thus we can write the Gaussian gap equation (2.31)in the neighbourhood of criticality in the form, 



(4.18) 



m w m 



I + 12A/ 3 (m w 0) - f^/l^rn^O)- -^-ln(l - e^'AV o }) + 



e 

7rL 2 



Mi-^ eL2 ^^) + — E 



7rLl 

e -eV / o(™~°)\/ L ?™?+ i 2™l 



"1 ,"2 = 1 



(4.19) 



Taking m 2 , = a(T/Tn — 1), with a > we obtain from the above equation for m 2 = 0, the critical temperature as a 
function of the wire transverse dimensions, L\ L 2 and the bulk transition temperature, T , 



T c (Li,L 2 ) — T Q 



+^_;„(l_ e =fav^nSSi))_HL y 

HIT J ,n /77T ^ 



anL 2 



an 



m ,n2=l 



(4.20) 



For reasons that will become apparent later, let us take (L 2 = pL\), L\ = L, with p a positive constant; then we 
get from Eqs.(4.13) and (4.15), 



(4.21) 



where 



27T V p ' 7T 



oo _^ 

^ (Ko(2irpnin2) + -K Q (2Trnin 2 /p)) 



_ni ,n2 = l 

and the critical temperature (4.20) is rewritten in the form, 



(4.22) 



T c (L,p) = T 



+ 



1 



12AO, 



+ 



" + J_j„(l_e- 



ejLC; 



aL 2na\fL airL 



p )+ 



airpL 



ln(l 



IhC. 



)-— y 



a7rL 



«1,"2 = 1 



-ey/LC p yJnl+p 2 



(4.23) 



This equation describes the behavior of the critical temperature of a type-I superconducting wire having a rect- 
angular cross-section, in terms of one of its transverse dimensions (L\ = L) and of the parameter p = LijL\. We 
recall that in this expression all gauge fluctuations have been taken into account by means of the Gaussian effective 
potential. These contributions are given by the last four terms in Eq,(4.23). 

Comparison to some experimental results: 



Up to now we have used Natural Units, c = K = k B = 1- To proceed, we restore SI units, remembering that 
a factor l/k B To is implicit in the exponent of equation (3.1). Next we rescale the fields and coordinates by <f> — ► 
4>new = \fio/k B To(j), A — ► A new = ^£,a/k B T Q A and x — ► x new = x/£o, with £ = 0A&hv F /k B T o , being the intrinsic 
coherent length and v F the Fermi velocity of a particular material. This implies that, in the Ginzburg-Landau model, 
the parameters a, A and e become dimensionless and are given by [14], 



a = l, A = 111.08(-^) 2 , e 

J-F 



2.59 



avp 



(4.24) 



where Tp and a are respectively the Fermi temperature and the fine structure constant. Thus, equation (4.23) is 
rewritten as, 



T c {L,pL) = T 



1 - 1333.0(|)(^) 2 Cj ; - 17.2.1:^:10 Vy'V ( Vw\ - 



+ 15.594 a T0- 3 (|) WFc ln(l e -°' 22128 V*V^ : ) + 



+ 15.594xlO- 3 (^)v Fc ln{l-e-°™ 
pL 



\J C p V Fc 



) + 



-0.031187(^)«f c J2 



0.22128y / ^ A /C p «F<:V™?+P 2 "l 



(4.25) 



As an application of this equation we plot the behavior of the critical temperature T c (L,p) as a function of 1/L to a 
superconducting wire made from niobium, characterized by vp = 1.37 x 10 6 m/s, To = 9.3 K, and Tp = 6.18 x 10 4 
K. We have considered three cases, with p = 1, p = 4 and p — 10 corresponding to wires of different rectangular 
cross-sections. 



10-, 




100/L(A D "' ) 



FIG. 1. Critical temperature T c as function of one of the transverse dimensions 1/L, from Eq.(4.25) for a superconducting 
wire made from Niobium of rectangular section cross A — L x pL. The figure shows the behavior of the critical temperature 
for the cases, p = 1, p = 4 and p — 10, respectively. 

In Fig.l, we show the quasi linear behavior of the critical temperature of a niobium wire as a function of L, for 
p = 1, p = 4 and p = 10. These curves suggest the existence of minimal values of L (or of the rectangular wire 
cross-sections), for which the surperconducting phase transition is suppressed. The approximate values for each case 
are: L™" « 2.5 nm, L™I% « 1.0 nm and Z/™" s=s 0.5 nm respectively. 
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FIG. 2. Critical temperature T c as function of the inverse thickness 1/L, for a Niobium film 

Also, it can be shown numerically that as we take larger and larger values of p (p — > oo) in Eq.(4.25) all the curves 
tend to coincide, corresponding to a sample in the form of a film of thickness L. A numerical analysis of Eq.(4.22) 
shows that as p — > oo, the quantity C p tends to a constant value, C p = Replacing this value in Eq.(4.25) we 
get the critical temperature for a film. In Fig. 2 we plot Eq.(4.25) in this limiting situation, which corresponds to a 
superconducting film of thickness L. In this case we find a minimal film thickness ss 0.13nm. 

V. CONCLUSIONS 

In this paper we have considered the confined Ginzburg-Landau model, in the transverse unitarity gauge, as a 
model to describe samples of superconducting materials not in bulk form. To generate the contributions from gauge 
fluctuations, we have used the Gaussian effective potential [2-5], which allows to obtain a gap equation that can be 
treated with the method of recent developments [8,9]. We have derived a critical equation that describes the changes 
in the critical temperature due to confinement. Contributions from the self interaction of the scalar field and from 
the gauge field fluctuations can be clearly identified. Our approach suggests a minimal size sustaining the existence 
of charged and non-charged transitions for both superconducting wires and films. Such a kind of result could be of 
practical interest to define limits of miniaturization of superconducting wires in manufacturing electronic circuits. It 
is interesting to compare the minimal allowed value of L w 2.5 nm for a wire having a square cross-section that we 
have obtained, to the experimentally observed superconducting state in nanowires having a cross-section diameter of 
w 10 nm [29] 

Also, the behavior described in Eq.(4.25) for p — > oo (a film) could be contrasted with the linear decreasing of the 
critical temperature with the inverse of the film thickness, that has been found experimentally in materials containing 
transition metals, for example, in PB [20], in W-Re alloys [21], in Nb [22,24,26], Mo-Ge [23]; for some of these cases, 
this behavior has been explained in terms of proximity, localization and Coulomb-interaction effects. We can clearly 
see from Eq.(4.25), as noticed above, that a linear decreasing of the critical temperature with the inverse of film 
thickness is recovered when we take e = 0. Even though, with e^Oa quasi-linear decreasing of T c with 1/L can be 
directly obtained from Eq.(4.25), since the terms coming from effects of gauge fluctuations are very small as compared 
to the term generated from the self coupling. Also a comparison may be done with recent theoretical results for type 
II superconductors [28], where a similar behavior of the critical temperature with the film thickness has been found 
for non-charged transitions. 

However, it should be noticed that we have used 3-dimcnsional values for the phcnomcnological parameters. This 
means that we should restrict ourselves only to relatively thin wires and films, which could be considered as essentially 
3-dimensional objects; very thin wires and films can not be physically accommodated in the context of our model. For 



instance, the value we have obtained for L^]^ (iV6) (w 0.13 nm) is of the order of magnitude of a few Bohr radii and 
should be considered to correspond to a 2-dimcnsional system. So, on physical grounds, it is beyond the domain of 
validity of the model. Nevertheless our results obtained from the "pure" GL model are in qualitative agreement with 
the experimentally observed behaviors mentioned above. It is worth to emphasize that the quasi-linear character of 
the decreasing of the transition temperature obtained in this paper, emerge solely as a topological effect of the spatial 
compactification of the Ginzburg-Landau model. 

Also, it should be remarked that our formalism can not account for microscopic details present in real samples of 
materials nor for the influence of manufacturing aspects, like the kind of the subtract on which a film is deposited, or 
the presence of impurities. This means that our expression for the transition temperature should be further improved 
(in particular the GL parameters), and more work is needed to better understand the behavior of real samples. This 
will be the subject of future investigation. 
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